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ABSTRACT 


This  Phase  I  effort  demonstrated  the  feasibility  of  using  scatterometry  to  provide  sensors 
for  process  control  and  overlay  measurement  for  structures  having  sub- pm  linewidths 
(CDs).  The  processes  which  we  investigated  were  etching  poly-Si  and  selective 
deposition  of  W,  with  particular  attention  directed  toward  application  in  a  cluster  tool 
environment.  We  analyzed  structures  having  CDs  of  ~0.20  pm,  ~0.35  pm,  and  ~0.50  pm. 
We  used  simulation  techniques  to  show  that  scatterometry  is  capable  of  monitoring  the 
width  and  height  of  the  lines  during  these  processes,  summarised  as  follows: 


Structure 

CD  Prediction  Error 
(me) 

Height  Prediction  Error 
(nm) 

131  »TH  3EEESE  J!EWi 

0.4 

8 

3 

21 

3 

14 

1 

4 

4 

5 

These  results  can  be  improved  by  implementing  the  sub-range  prediction  techniques 
described  We  have  demonstrated  the  ability  of  scatterometry  to  provide  overlay 
measurement  in  a  simple  manner  which  is  summarized  as  follows: 

OVERLAY  MEASUREMENT  ERROR:  <;  0.5  nm 

Note  that  we  have  applied  these  same  simulation  techniques  in  similar  applications 
involving  CDs  in  the  range  between  0  1 2  pm  and  2  pm,  and  results  have  agreed  well  with 
experimental  results.  In  a  '*hasc  JI  effort  wc  would  verify  these  simulations  by 
constructing  scatterometers  and  implementing  them  in  reactors. 
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1.0  INTRODUCTION 


The  development  of  critical  dimension  (CD)  metrology  techniques  has  become  crucial 
as  the  CDs  of  semiconductor  devices  have  shrunk  toward  0.25|im.  The  tolerances  of 
the  fabrication  process,  or  error  budget,  are  approximately  10%  of  the  CD.  The 
shrinking  tolerances  have  led  to  the  industry’s  steady  drift  from  optical  microscopy  to 
scanning  electron  microscopy  (SEM)  |1|  as  the  primary  reference  method  for  CD 
metrology.  Additional  CD  metrology  tools  include  scanning  force  microscopy  (SFM) 
[2]  and  electrical  linewidth  measurements  [3].  The  critical  performance  issues  with 
every  metrology  instrument  are  calibration,  accuracy  and  precision  |4|.  C  her  desirable 
attributes  are  that  the  tool  be  non-contact,  non-destructive,  rapid  and  applicable  itt-si tu. 
However,  these  three  metrology  techniques  suffer  from  some  combination  of  being 
slow,  destructive,  tedious,  and  unsuitable  for  in-situ  application. 

Optical  scatterometry  applied  to  periodic  structures  has  all  of  the  above  attributes,  and 
hence  has  tremendous  potential  as  a  CD  metrology  tool.  Scatterometry  consists  of  the 
measurement  and  analysis  of  light  scattered  from  a  sample  surface  J5].  In  the  case  of  a 
periodic  structure,  such  as  arrays  or  repetitive  cells,  the  light  is  diffracted  orders  at 
angular  locations  specified  by  the  grating  equation 

sin@.  +  sin©n  =  —  ,  (1) 

d 

where  0/  is  the  angle  of  incidence,  in  the  angular  location  of  the  n-th  diffraction 
order,  n  is  an  integer,  X  is  the  wavelength  of  incident  radiation,  and  d  is  the  period 
(pitch)  of  the  structure.  Figure  1  illustrates  the  geometry  used  in  the  description  of 
diffraction. 


Figure  1 :  Diffraction  geometry  used  in  this  effort 
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Information  about  the  diffracting  structure  is  obtained  by  interpreting  the  distribution  of 
intensities  in  the  various  diffracted  orders.  This  distribution  is  extremely  sensitive  to  the 
width,  height,  and  sidewall  profile  of  the  diffracting  structure.  Scatterometry  is  an 
optical  metrology  technique  that  does  not  suffer  from  the  limitations  of  optical 
microscopy  because,  unlike  imaging,  the  diffracted  light  is  not  combined  or  imaged 
prior  to  analysis. 

When  measuring  diffraction  from  gratings,  the  polarization  of  the  incident  beam 
measurably  affects  the  intensity  of  the  diffracted  beams.  For  TE  polarization,  the  incident 
F,  field  is  parallel  to  the  grating  lines  When  the  E  field  is  perpendicular  to  the  grating 
lines,  the  beam  is  TM  polarized.  Typically  we  use  TE  polarization,  although  the 
polarization  may  be  changed  in  a  controlled  manner  to  increase  the  amount  of  information 
obtained  from  a  measurement 

For  a  fixed  angle  of  incidence  of  the  laser  beam,  the  number  of  diffracted  orders  is 
determined  by  the  ratio  A./d.  The  larger  this  ratio  is,  the  greater  the  number  of 
diffraction  orders,  and  consequently,  the  greater  the  information  content  of  the 
scatterometry  data. 

For  small  pitch  structures  in  which  the  period  of  the  grating  is  comparable  to  the 
wavelength  of  the  incident  probe  beam  as  illustrated  in  Figure  1,  only  a  few 
propagating  diffraction  orders  exist.  When  the  angle  of  incidence  is  fixed,  little 
information  is  available  from  a  scatterometer  measurement,  and  alternative 
measurement  arrangements  must  be  implemented.  Either  the  angle  of  incidence  or  the 
wavelength  of  the  incident  beam  (or  both)  can  be  scanned,  and  the  intensity  of  one  or 
more  diffraction  orders  can  be  continuously  monitored.  For  this  Phase  I  effort  we 
demonstrated  the  feasibility  of  scatterometers  which  incorporate  these  two  techniques; 
this  is  discussed  in  Section  2. 

The  behavior  of  the  diffracted  fields  and  their  dependence  on  the  diffracting  structure 
lineshape  can  be  understood  using  a  rigorous  form  ism.  A  formalism  we  have  used 
with  great  success  to  simulate  and  interpret  o  i  fraction  data  encountered  using 
scatterometry  is  rigorous  coupled  wave  theory  (RCWT)  |6|.  RCWT  numerically  solves 
Maxwell’s  equations  for  the  case  of  an  incident  wave  illuminating  a  periodic  structure. 
It  is  a  first-principles  solution  of  the  diffraction  problem  and  is  therefore  quantitative 
and  accurate.  Use  of  RCWT  allows  for  the  theoretical  solution  of  the  forward  grating 
diffraction  problem,  that  of  predicting  the  distribution  of  intensities  in  the  various 
diffraction  orders  when  a  laser  beam  is  incident  on  a  grating  of  known  shape. 

The  inverse  grating  diffraction  problem ,  predicting  the  lineshape  of  the  diffracting 
structure  from  knowledge  of  the  diffracted  light  characteristics,  is  quite  difficult.  This 
is  because  the  diffraction  characteristics  vary  in  a  complex,  non-intuitive  manner  as  the 
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lineshape  of  the  diffracting  structure  is  changed.  The  inverse  problem  has  been  a  topic 
of  active  research  for  many  years.  Our  approach  to  this  has  involved  multivariate 
statistical  techniques  (e  g.  Chemometrics)  and  has  provided  very  good  results.  These 
same  statistical  techniques  have  been  applied  in  this  Phase  1  effort  and  are  described  in 
Section  2.  Diffracted  Scatter  Analysis  (DSA)  (7),  used  for  the  solution  of  the  inverse 
problem,  combines  scatterometry,  theoretical  simulation,  and  statistical  analysis  to 
predict  critical  dimensions  from  a  measured  optical  scatter  pattern. 
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2.0  OVERVIEW  OF  OUR  APPLICATION  OF  SCATTEROMETRY  TO 
PROCESS  CONTROL  AND  OVERLAY  MEASUREMENT 


As  explained  previously,  for  a  diffracting  structure  in  which  the  pitch  (and  CD)  of  the 
diffracting  structure  is  small  compared  to  the  incident  probe  beam  wavelength,  the 
angle  of  incidence  or  the  wavelength  can  be  varied  in  order  to  obtain  a  large  amount  of 
information  in  a  measurement.  A  "2-0"  scatterometer  is  the  name  applied  to  the 
arrangement  in  which  the  angle  of  incidence  is  scanned,  and  a  "X-scan"  scatterometer  is 
the  name  applied  to  the  other  situation  in  which  the  wavelength  of  the  incident  probe 
beam  is  scanned.  In  these  two  situations,  the  shapes  of  the  curves  that  are  obtained  are 
extremely  sensitive  to  the  lineshape  of  the  diffracting  structure.  Figures  2  and  3 
illustrate  the  two  scatterometer  arrangements.  In  a  Phase  II  effort  we  would  construct 
and  implement  practical  versions  of  these  scatterometers  which  have  no  moving  parts 
and  other  improvements  which  make  them  more  amenable  to  fab  line  (cluster  tool)  use. 

CENTER  OF 


Figure  2  A  2-0  scatterometer  which  measures  the  power  of  a  diffracted  order  as  a 
function  of  incident  angle 


Figure  3  A  X-scan  scatterometer  which  measures  the  power  of  a  diffracted  order  as  a 
function  of  incident  wavelength. 
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We  first  identified  two  critical  processes  used  in  a  cluster  tool  environment  that  require 
sensors  for  process  control:  poly-Si  etch  and  selective  deposition  of  W;  additional 
details  for  this  decision  are  discussed  below.  The  RCWT  approach  was  used  to 
theoretically  generate  the  scatter  patterns  corresponding  to  the  2-0  and  X-scan 
scatterometer  arrangements  over  a  parameter  space  defined  by  expected  variation  of 
linewidth  and  height  of  the  structure.  The  wavelength  used  for  the  probe  beam  for  the 
0.20  pm  geometry  was  325  nm,  and  633  nm  was  used  for  the  0.35  pm  and  0.50  pm 
geometries.  The  effects  of  process  variations  were  simulated  by  allowing  the  height 
and  width  of  the  lines  to  vary  around  their  nominal  values. 

A  similar  approach  was  taken  for  overlay  measurement.  The  structure  which  was 
investigated  consisted  of  two  gratings.  One  grating  had  Si  lines  and  Si02  spaces,  and 
the  other  grating  consisted  of  developed  photoresist  and  was  located  on  top  of  the 
Si/SiO?  grating.  The  two  gratings  were  modeled  to  have  the  same  pitch  but  were 
offset  from  each  other  by  different  amounts.  The  wavelength  of  the  probe  beam  used 
in  these  simulations  was  633  nm.  Additional  details  are  given  in  Section  4. 

The  theoretical  scatter  patterns  provided  the  calibration,  or  reference  set  of  data  needed 
to  train  the  multivariate  statistical  analysis  technique.  This  was  used  to  solve  the 
inverse  problem  of  predicting  the  lineshape  and  overlay.  We  used  partial  least  squares 
(PLS)  for  this  purpose  |7,8|. 

Cross  validation  was  used  to  quantify  the  expected  performance  of  a  scatterometer  and 
analysis  technique  when  applied  to  a  specific  structure  Cross  validation  is  the  process  of 
repeatedly  calibrating  the  statistical  technique  using  a  subset  of  the  calibration  data,  using 
the  rest  of  the  calibration  data  as  "unknowns"  which  are  submitted  to  the  statistical 
technique  for  prediction,  and  recording  the  prediction  error  for  each  "unknown"  The 
standard  deviation  of  the  prediction  errors  is  the  standard  error  of  prediction  (SEP)  of  the 
proposed  instrument;  i.e.  the  SEP  values  listed  in  the  sections  below  represent  the  l-o 
precision  of  predicting  the  parameters  of  interest 

We  calculated  the  SEP  for  three  different  noise  situations  for  most  applications.  The  first 
situation  assumes  that  both  the  calibration  data  and  the  measurement  data  are  noise-free. 
The  second  situation  assumes  that  only  the  measurement  data  contains  noise.  The  third 
situation  assumes  that  both  the  measurement  data  and  the  calibration  data  contain  noise. 
Examining  these  three  cases  characterizes  the  noise  immunity  of  the  proposed  instrument. 

Our  experience  over  the  past  five  years  investigating  the  diffraction  from  etched  structures 
and  other  relief  gratings  has  been  that  the  0-order  and  Ist-order  diffracted  beams  are 
typically  strong  and  able  to  be  seen  with  the  unaided  eye.  Because  of  this  the 
measurements  contain  low  noise  levels,  and  a  conservative  simulation  involves  the  second 
case  described  above.  The  simulated  calibration  data  is  practically  noise  free. 


8 


The  noise  model  chosen  was 


I  =  I0  +  k*N*sqrt(I0)  ,  (2) 

where  N  is  a  zero  mean  unit  variance  Gaussian  random  variable,  k  is  a  constant  modifier 
used  to  change  the  noise  power,  Io  is  the  noise  free  measurement,  and  I  is  the  noisy 
measurement.  This  noise  model  mimics  electronic  noise  in  the  measurement  system.  The 
value  k=0.002  was  chosen  to  heavily  degrade  the  weak  signals  (with  Iq  ~  0.000008)  in  the 
simulated  measurement  while  leaving  the  stronger  signals  relatively  intact.  This  added 
simulated  dynamic  range  limitations  to  the  experiments. 

This  app/oach,  i.e.  diffraction  simulation  using  the  RCWT  coupled  with  multivariate 
statistical  analysis  (PLS),  provides  an  excellent  platform  for  instrument  design.  Simulation 
can  precisely  characterize  the  behavior  of  a  prototype  instrument  before  the  instrument  is 
fabricated.  The  characteristics  of  a  proposed  instrument,  including  the  measurement 
accuracy  and  precision,  are  known  prior  to  construction  An  additional  benefit  of  this 
method  of  design  is  that  the  instrument  is  calibrated  in  the  design  phase.  Because  the 
simulation  program  is  based  on  first  order  principals,  the  final  instrument  calibration  is  well 
characterized  and  is  independent  of  any  commercially  available  calibration  standard.  We 
have  used  this  same  approach  in  similar  applications  in  which  the  calculations  and 
experimental  verification  were  in  close  agreement,  such  as  characterizing  the  lineshape  of 
developed  photoresist.  See  Reference  5,  which  is  attached  in  Appendix  1. 


9 


3.0  CRITICAL  DIMENSION  METROLOGY 


The  first  task  we  performed  was  to  determine  relevant  applications  for  scatterometry  in 
today's  and  tomorrow's  microelectronics  industry.  We  visited  Applied  Materials  in  Santa 
Clara,  CA,  TI  in  Dallas,  TX,  SC  Technologies  in  Livermore,  CA,  and  Sandia  National 
Laboratories  in  Albuquerque,  NM,.  Representatives  of  Applied  Materials  and  Tencor 
subsequently  visited  us  in  Albuquerque.  We  communicated  with  Novellus  in  Santa  Clara, 
CA,  Drytech  in  Wilmington,  MA,  and  AT&T  in  Allentown,  PA.  These  contacts  reported 
or  confirmed  that  etching  poly-Si  and  filling  trenches  with  W  are  presently  relevant 
processes  that  would  benefit  from  advanced  metrology. 


Figure  4.  The  five  structures  that  were  examined  to  demonstrate  the  applicability 
of  scatterometry  to  monitor  line  shapes  of  etched  poly-Si  (parts  a,  b,  and  c)  and 
deposited  W  (parts  d  and  e). 


Figure  4  shows  the  five  specific  structures  we  examined  which  correspond  to  two 
processes,  with  structures  having  nominal  CDs  of  -0.20  pm,  -0  35  pm,  and  -0.50  pm. 
The  poly-Si  etching  involved  a  photoresist  mask  that  was  0.5  pm  thick  for  the  0  2  pm  CD 
structures,  and  1.0  pm  thick  for  the  0.35  pm  and  0.50  pm  structures.  The  selective  W 
deposition  process  involved  filling  trenches  in  a  SiC>2  grating  over  a  Si  substrate.  The 
structures  shown  in  Figures  4d  and  4e  show  the  trench  partially  filled  with  W.  The  second 
process  was  etching  poly-Si  through  a  photoresist  mask. 


The  line  width  of  the  structures  was  varied  by  ±1 0%  of  the  nominal  value  to  accommodate 
for  process  variations  that  might  occur.  In  addition,  for  each  structure  the  size  of  the  lines 
and  spaces  was  approximately  the  same.  This  is  realistic  if  issues  involving  proximity 
effects  are  to  be  avoided. 

We  emphasized  0.35  pm  and  0.50  pm  geometries  in  this  effort  because  of  their  relevance 
to  our  Phase  II  effort  which  we  will  propose.  Specifically,  as  discussed  below,  we  will 
propose  to  implement  a  scatterometer  on  a  deposition  or  etch  system  to  monitor  the 
process,  and  0.20  pm  geometries  are  not  compatible  with  the  systems  available  to  us  In 
addition,  the  larger  geometries  represent  a  logical  next  step  in  the  progression  from 
current  technology  to  smaller  geometries.  We  will  approach  measuring  0.20  pm 
structures  ex-situ  as  discussed  in  Section  7. 

Figures  5,  6,  and  7  show  examples  of  data  which  we  generated  using  RCWT.  Figure  5 
corresponds  to  the  structure  of  Figure  4a  (etched  poly-Si  and  CD  ~  0.20  pm)  and  shows 
the  variation  of  the  lst-order  diffracted  power  with  incident  angle  of  the  probe  beam  for 
three  different  line  widths.  Figure  6  corresponds  to  the  structure  of  Figure  4b  (etched 
poly-Si  and  CD  ~  0.35  pm)  and  shows  the  variation  of  the  1  st-order  diffracted  power  with 
incident  wavelength  of  the  probe  beam  for  three  line  widths.  Fij,  ire  7  corresponds  to  the 
structure  of  Figure  4d  (selective  W  deposition  and  CD  ~  0.35  pm)  and  shows  the  variation 
of  0-order  power  with  incident  angle  of  the  probe  beam  for  three  line  heights  The 
differences  between  the  three  curves  in  each  of  the  figures  illustrates  the  sensitivity  of  the 
scatterometry  technique  to  parameter  changes,  and  this  is  what  partial  least  squares 
exploits  for  prediction.  Note  that  only  three  curves  are  included  in  each  figure  to  avoid 
confusion  that  would  occur  if  all  the  data  were  plotted. 


JO  *o 

I  nc  1  d  •n  t  Ang  I  • 


LEGEND 

M-W-M  0.  IB  0  O  0  0.19  »  *•  X  0  20 

Figure  5:  Simulated  2-®  scatterometer  data  showing  the  1  st-order  diffraction 
from  characterizing  an  etched  poly-Si  structure  of  CD  -  0.2  pm,  for  three  different 
line  widths;  ledgend  units  are  pm 
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Figure  6:  Simulated  X-Scan  scatterometer  data  showing  the  Ist-order  diffraction 
from  characterizing  an  etched  poly-Si  structure  of  CD  ~  0.3S  pm,  for  three 
different  line  widths;  ledgend  units  are  pm. 
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Figure  7:  Simulated  2-0  scatterometer  data  showing  0-order  diffraction  from 

■  characterizing  selective  W  deposition  structure  of  CD  -  0.35  pm,  for  three 

different  line  heights;  ledgt  nd  units  are  pm. 
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3.1  Si  Etch,  0.20  pm  CD 


The  0.20/0.4  poly-Si  etch  structure  is  shown  in  Figure  4a  which  consists  of  a  0.5  pm  thick 
photoresist  grating  mask  over  the  poly-Si  layer.  The  wavelength  of  the  incident  probe 
beam  was  325  nm  and  TE  polarization  was  used.  The  simulation  started  with  a  line  height 
of  0.00  pm  and  continued  until  a  line  height  of  0.40  pm  was  reached,  increasing  in 
increments  of  0.02  pm,  to  simulate  a  total  of  21  heights.  The  line  width  was  varied  by 
±10%  around  the  nominal  width  of  0.20  pm,  in  steps  of  0.01  pm  for  a  total  of  5  widths. 
For  this  geometry  the  total  number  of  simulations  was  21x5=105.  Each  simulation 
consisted  of  a  2-0  scan  over  the  range  of  0°  io  78°  for  the  angle  of  incidence,  taken  in 
steps  of  2°.  In  this  application  we  included  noise  in  the  measurement  data  and  the 
calibration  data  (i.e.,  the  third  noise  situation  discussed  in  Section  2). 

Table  1  shows  the  performance  of  the  2-0  technique  is  very  good  at  predicting  the  line 
width  and  the  height  To  improve  on  the  height  prediction,  we  implemented  a  technique 
called  "sub-range  analysis".  This  strategy  first  utilizes  the  full-range  results  of  Table  1  to 
locate  the  line  height  to  within  14  nm.  This  places  the  height  into  one  of  three  sub-ranges: 
0  nm  to  200  nm,  100  nm  to  300  nm,  and  200  nm  to  400  nm  The  predictions  for  line 
height  are  then  performed  over  the  smaller  sub- range  to  provide  a  significant  improvement 
in  the  height  SEP,  as  illustrated  in  Table  2. 


Measurement  Technique:  2-0 

Diffraction  Order(s)  Used: 

0&1 

Width  SEP  (nm) 
noise  in  measurement 

0.4 

Height  SEP  (nm) 
noise  in  measurement 

14 

Table  1 .  The  SEPs  for  measuring  critical  dimensions  of  the  0.20/0.4  etched  poiy- 
Si  structure  using  the  full  range  of  line  height  for  prediction 


Sub-Range  (nm): 

0- 

200 

100  - 
30) 

200- 

400 

1.3 

2.5 

8.5 

Table  2  Sub-range  analysis  SEPs  for  line  height  prediction  using  unknown  line 
width  of  the  0.20/0  40  etched  poly-Si  structure  A  2-0  measurement  technique 
and  both  diffraction  orders  were  used  for  the  analysis. 

Note  that  this  refinement  in  the  analysis  process  does  not  involve  any  additional 
measurements  and,  because  all  the  calibrations  are  already  performed,  can  be  computed 
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quickly.  After  training  the  PLS  technique,  the  total  analysis  time  required  for  prediction  is 
on  the  order  of  a  few  msec.  The  sub-range  process  can  be  repeated  to  further  improve 
prediction  results.  We  are  very  pleased  with  this  approach  and  have  applied  it  to  some  of 
the  other  structures  discussed  below. 


3.2  Si  Etch,  0.35  pm  CD 

The  0.35/0.7  poly-Si  etch  structure  is  shown  in  Figure  4b.  The  process  simulation  began 
with  a  1 .0  pm  thick  photoresist  grating  over  a  poly-Si  layer.  The  end  point  was  when  the 
poly-Si  had  been  etched  0  5  pm  deep.  Simulations  were  run  starting  at  the  beginning  and 
at  .02  pm  intervals  until  the  0.5  pm  line  height  was  obtained.  The  total  number  of  etch 
heights  was  26.  The  line  width  was  varied  ±10%  around  the  nominal  width.  For  the 
0.35/0.7  structure,  the  range  was  0  315  pm  to  0.385  pm.  A  step  size  of  0.007  pm 
resulted  in  1 1  different  line  widths.  The  total  number  of  simulations  was  26x11=281.  The 
wavelength  of  the  incident  probe  beam  was  633  nm 


Measurement  Technique: 

2-0 

X-scan 

Diffraction  Order(s)  Used: 

0-order 

1 -order 

0&1 

0-order 

1 -order 

0&1 

Width  SEP  (nm) 
no  noise 

2 

5 

1 

10 

11 

3 

Width  SEP  (nm) 
noise  in  measurement 

2 

5 

1 

10 

12 

3 

Width  SEP  (nm) 
noise  in  measurement  & 
noise  in  calibration 

3 

6 

1 

10 

11 

3 

Height  SEP  (nm) 
no  noise 

42 

47 

29 

51 

27 

20 

Height  SEP  (nm) 
noise  in  measurement 

44 

57 

29 

53 

28 

21 

Height  SEP  (nm) 
noise  in  measurement  & 
noise  in  calibration  (nm) 

43 

56 

on 

<C7 

52 

27 

Table  3.  The  SEPs  (standard  errors  of  prediction)  for  measuring  critical 
dimensions  of  the  0.35/0  70  etched  poly-Si  structure. 

The  results  shown  in  Table  3  indicate  that  the  0.35/0.7  Si  etch  structure  can  be  predicted 
accurately  using  either  arrangement,  although  the  X-scan  technique  provides  slightly  better 
results  The  slight  difference  in  the  SEP  of  the  predicted  width  from  the  two  arrangements 
is  inconsequential  The  performance  of  both  techniques  could  be  improved  using  the  sub¬ 
range  technique  described  above  in  connection  with  the  0.20  pm  structure. 
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3.3  Si  Etch,  0.50  pm  CD 

The  0.5/1. 0  poly-Si  etch  structure  is  shown  in  Figure  4c.  The  process  simulation  began 
with  a  1 .0  pm  thick  photoresist  grating  over  a  layer  cf  poly-Si,  representing  a  0.0  pm  deep 
etch.  Simulations  were  run  at  .02  pm  intervals  until  a  line  height  of  1 .0  pm  was  reached. 
The  total  number  of  line  heights  was  51.  The  line  width  was  varied  ±10%  around  the 
nominal  width.  For  the  0.5/1. 0  ructure,  the  range  was  0.45  pm  to  0.55  pm.  A  step  size 
of  0.01  pm  resulted  in  1 1  different  line  widths.  The  total  number  of  simulations  was 
51x11=561.  The  wavelength  of  the  incident  probe  beam  was  63  3  nm. 


Measurement  Technique: 

2-0 

A.-scan 

Diffraction  Orderfs)  Used: 

0-order 

1 -order 

0&1 

0-order 

1  -order 

0&1 

Width  SEP  (nm) 
no  noise 

5 

5 

2 

6 

8 

3 

Width  SEP  (nm) 
noise  in  measurement 

5 

6 

2 

6 

8 

3 

Width  SEP  (nm) 
noise  in  measurement  & 
noise  in  calibration 

5 

5 

2 

6 

8 

3 

Height  SEP  (nm) 
no  noise 

241 

257 

209 

278 

273 

227 

Height  SEP  (nm) 
noise  in  measurement 

244 

277 

215 

279 

273 

227 

Height  SEP  (nm) 
noise  in  measurement  & 
noise  in  calibration  (nm) 

242 

273 

215 

279 

273 

229 

Table  4.  The  SEPs  for  measuring  critical  dimensions  of  the  0.5/1 .0  poly-Si  etch 
structure. 


As  Table  4  shows,  the  line  width  prediction  results  were  excellent,  but  the  height  estimates 
were  poor  in  ail  cases.  Examination  of  the  data  revealed  a  periodic  variation  of  the 
diffracted  intensity  with  increasing  line  height.  The  variation  is  an  interferometric  effect  of 
the  optical  path  difference  inside  and  outside  the  trench.  The  PLS  analysis  performs  best 
when  the  measurement  varies  linearly  with  the  critical  dimension,  and  the  solution  quickly 
degrades  when  this  is  not  the  case  The  PLS  analysis  is  also  poor  when  the  measurement 
is  not  unique  over  the  range  of  critical  dimensions  under  study,  as  occurs  here  with 
periodic  variations. 

We  then  applied  sub-range  analysis  described  in  Section  3.1  to  improve  the  line  height 
predictions.  The  height  was  divided  into  5  overlapping  sub-ranges  between  0.0  pm  and 
1.0  pm,  and  the  width  was  divided  into  4  overlapping  sub-ranges  between  0.45  pm  and 
0  55  urn.  The  full-range  results  summarized  in  Table  4  were  used  to  locate  the  prediction 


into  one  of  the  20  "bins"  corresponding  to  a  particular  combination  of  width  and  height 
sub-ranges.  We  then  applied  PLS  to  determine  the  SEP  for  each  sub-range.  We  used  0- 
order  and  lst-order  data  in  which  the  measurement  data  and  the  calibration  data  had  noise 
(i.e.  the  third  noise  situation). 

Tables  5  and  6  summarize  the  line  height  SEP  obtained  for  each  of  the  20  "bins"  using  the 
2-0  technique  and  the  A.- scan  technique,  respectively.  It  can  be  seen  that  the  height 
prediction  was  significantly  improved  to  provide  SEPs  that  are  less  than  14  nm  for  the  X- 
scan  technique  and  less  than  36  nm  for  the  2-0  technique  Thus  the  A.-scan  technique  is 
considered  to  provide  the  better  overall  performance  for  this  specific  application.  We  are 
very  pleased  with  both  of  these  results. 


Linewidth  & 
Height  Ranges 

0.45  -  0.49 
(pm) 

0.47-0.51 

(pm) 

0.49-0.53 

(pm) 

0.51-0.55 

(pm) 

0.0-0.36 

(pm) 

0.0053  pm 

0.0061  pm 

0.0083  pm 

0.0086  pm 

0.18-0.54 

(pm, 

0.0185  pm 

0.0187  pm 

0.0280  pm 

0.0238  pm 

0.36-0.72 

(pm) 

0.0220  pm 

0.0255  pm 

0.0225  pm 

0.0213  pm 

0.54  -  0.90 
(pm) 

0.0346  pm 

0.0276  pm 

0.0312  pm 

0.0362  pm 

0.72-  1.00 
(pm) 

0.0199  pm 

0.0061  pm 

0.0138  pm 

0.0099  pm 

Table  5.  Sub-range  analysis  SEPs  for  line  height  prediction  using  unknown  line 
width  of  the  0.5/1  0  etched  poly-Si  structure  n  2-0  measurement  technique 
and  both  diffraction  orders  were  used  for  the  analysis. 


Linewidth  & 
Height  Ranges 

0.45  -  0.49 
(pm) 

0.47-0.51 

(pm) 

0.49-0.53 

(pm) 

0.51-0.55 

(pm) 

0.0-0.36 

(pm) 

0.0037  pm 

0.0043  pm 

0.0061  pm 

0.0088  pm 

0.18-0.54 

(pm) 

0.0094  pm 

0.0086  pm 

0.0137  pm 

0.0118  pm 

0.36  -  0.72 
(pm) 

0.0100  pm 

0.  X)83  pm 

0.0143  pm 

0.0115  pm 

0.54  -  0.90 
(pm) 

0.0145  pm 

0.0135  pm 

0.0101  pm 

0.0095  pm 

0.72-1.00 

(pm) 

0.0036  pm 

- - - 

0.0031  pm 

0.0037  pm 

0.0036  pm 

Table  6.  Sub-range  analysis  SFPs  for  line  height  prediction  using  unknown  line 
width  of  the  0.5/1 .0  etched  poiy-Si  structure.  A  X-scan  measurement  technique 
and  both  diffraction  orders  were  used  for  the  analysis. 


3.4  Selective  W  Deposition,  0.35  pm  CD 

The  0.35/0.7  W  deposition  structure  is  shown  in  Figure  4d.  The  process  simulation  began 
with  a  0.7  pm  pitch,  0.5  pm  deep  Si(>2  grating  over  a  bare  Si  wafer,  representing  a  0  0 
pm  W  fill.  Simulations  were  run  starting  at  0.0  pm  and  at  .02  pm  intervals  until  the  0.5 
pm  deep  trench  was  completely  filled.  The  total  number  of  deposition  heights  was  26. 
The  line  width  was  also  varied  ±10  %  around  the  nominal  width  to  provide  a  range  of 
0.315  pm  to  0.385  pm,  A  step  size  of  0.007  pm  resulted  in  )  I  different  line  widths.  The 
total  number  of  simulations  was  26  x  11=286.  The  wavelength  of  the  incident  probe  beam 
was  633  run 
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Measurement  Technique: 

2-0 

X-scan 

Diffraction  Order(s)  Used: 

0-order 

1 -order 

0&1 

0-order 

1 -order 

0&1 

Width  SEP  (nm) 
no  noise 

7 

6 

2 

1 

2 

<1 

Width  SEP  (nm) 
noise  in  measurement 

8 

7 

2 

2 

3 

1 

Width  SEP  (nm) 
noise  in  measurement  & 
noise  in  calibration 

8 

6 

2 

2 

3 

1 

Height  SEP  (nm) 
no  noise 

6 

16 

4 

5 

6 

2 

Height  SEP  (nm) 
noise  in  measurement 

8 

24 

4 

8 

15 

4 

Height  SEP  (nm) 
noise  in  measurement  & 
noise  in  calibration  (nm) 

8 

24 

4 

8 

16 

3 

Table  7.  The  jEPs  (standard  errors  of  prediction)  for  measuring  critical 
dimensions  of  the  0.35/0.7  W  deposition  structure. 

Table  7  illustrates  the  predicted  results  for  the  0.35/0  7  W  deposition  structure  are 
exceptional  for  both  scatterometer  arrangements.  The  X-scan  technique  slightly 
outperformed  the  2-0  technique  for  both  the  line  'vidth  and  height  predictions  However, 
these  differences  are  considered  to  be  inconsequential.  The  0* order  column  for  the  X-scan 
arrangement  suggests  using  a  single  fixed  detector  for  the  measurement,  and  this 
represents  a  simple  instrument  with  no  moving  parts  and  easy  alignment 

3.5  Selective  W  Deposition,  0.50  pm  CD 

The  0. 5/1.0  W  deposition  structure  is  shown  in  Figure  4e  The  process  began  with  a  1.0 
pm  pitch,  1.0  pm  deep  SiC>2  grating  over  a  bare  Si  wafer,  representing  a  0  0  pm  W  fill 
Simulations  were  run  at  0  0  pm  and  at  .02  pm  intervals  until  the  1.0  pm  deep  trench  was 
completely  filled.  The  total  number  of  deposition  heights  was  51  The  fine  width  was 
varied  ±10  %  around  the  nominal  width.  For  the  0.5/1 .0  structure,  the  range  was  0.45  pm 
to  0  55  pm  A  step  size  of  0.01  pm  resulted  in  1 1  different  line  widths.  The  total  number 
of  simulations  was  51x1 1=561.  The  wavelength  of  the  incident  probe  beam  was  633  nm. 

Table  8  illustrates  how  the  line  width  of  the  0. 5/1.0  W  deposition  structure  can  be  quite 
well  measured  using  either  technique,  although  the  X-scan  technique  performs  slightly 
better.  The  deposition  height  SEP  for  the  2-0  arrangement  is  marginal  because  the 
measurement  SEP  is  comparable  to  the  precision  available  using  SEM  characterization. 
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Measurement  Technique: 

2-0 

X-scan 

Diffraction  Order(s)  CJsed: 

0-order 

1 -order 

0&1 

0-order 

1 -order 

0&1  i 

Width  SEP  (nm) 
no  noise 

16 

25 

9 

12 

11 

4 

Width  SEP  (nm) 
noise  in  measurement 

16 

25 

9 

12 

13 

4 

Width  SEP  (nm) 
noise  in  measurement  & 
noise  in  calibration 

16 

25 

9 

12 

13 

5 

Height  SEP  (nm) 
no  noise 

67 

116 

49 

134 

125 

114 

Height  SEP  (nm) 
noise  in  measurement 

67 

124 

48 

134 

125 

114 

Height  SEP  (nm) 
noise  in  measurement  & 
noise  in  calibration  (nm) 

67 

148 

49 

134 

125 

114 

Table  8.  The  SEPs  (standard  errors  of  prediction)  for  measuring  critical 
dimensions  of  the  0.5/1 .0  W  deposition  structure 


We  applied  the  sub-range  analysis  described  in  Section  3  1  to  improve  the  line  height 
predictions.  The  height  was  divided  into  5  overlapping  sub-ranges  between  0.0  pm  and 
1.0  pm,  and  the  width  was  divided  into  4  overlapping  sub-ranges  between  0  45  pm  and 
0  55  pm  The  full-range  results  summarized  in  Table  8  were  used  to  locate  the  prediction 
into  one  of  the  20  "bins"  corresponding  to  a  particular  combination  of  width  and  height 
sub-ranges.  We  then  applied  PLS  to  determine  the  SEP  for  each  sub-range  We  used  0- 
order  and  1  st-order  data  in  which  the  measurement  data  and  the  calibration  data  had  noise 
(i  e.  the  third  noise  situation) 

Tables  0  and  10  summarize  the  line  height  SEP  obtained  for  each  of  the  20  "bins"  using 
the  2-0  technique  and  the  k-scan  technique,  respectively.  It  can  be  seen  that  the  height 
prediction  was  significantly  impi  oved  to  provide  SEPs  that  are  less  than  5  nm  for  the  X- 
3can  technique  and  less  than  6  nm  for  the  2-0  technique.  We  are  very  pleased  with  both 
of  these  results  and  consider  both  techniques  capable  of  characterizing  the  0.5/1 .0  W 
structure.. 
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Linewidtn  & 
Height  Ranges 

0.45-0.49 

(pm) 

0.47-0.51 

(pm) 

0.49-0.53 

(pm) 

0.51-0.55 

(pm) 

0.0-0.36 

(pm) 

0.0025  pm 

0.0027  pm 

0.0028  pm 

0.0017  pm 

0.18-0.54 

(pm) 

0.0040  pm 

0.0051  pm 

0.0047  pm 

0.0045  pm 

0.36-0.72 

(pm) 

0.0042  pm 

0.0034  pm 

0.0032  pm 

0.0040  pm 

0.54  -  0.90 
(pm) 

0.0026  pm 

0.0025  pm 

0.0025  pm 

0.0031  pm 

0.72  -  1.00 
(pm) 

0.0024  pm 

0.0018  pm 

0.0020  pm 

0.0021  pm 

Table  9.  Sub-range  analysis  SEPs  for  line  height  prediction  using  unknown  line 
width  of  the  0.5/1 .0  W  deposition  structure.  A  2-0  measurement  technique 
and  both  diffraction  orders  were  used  for  the  analysis. 


Linewidth  & 
Height  Ranges 

0.45-0.49 

(pm) 

0.47-0.51 

(pm) 

0.49-0.53 

(pm) 

0.51-0.55 

(pm) 

0.0-0.36 

(pm) 

0.0030  pm 

0.0056  pm 

0.0058  pm 

0.0038  pm 

0.18-0.54 

(pm) 

0.0030  pm 

0.0022  pm 

0.0019  pm 

0.0030  pm 

0.36-0.72 

(pm) 

0.0032  pm 

0.0019  pm 

0.0029  pm 

0.0031  pm 

0.54  -  0.90 
(pm) 

0.0028  pm 

0.0028  pm 

0.0026  pm 

0.0028  pm 

0.72-1.00 

(pm) 

0.0016  pm 

0.0019  pm 

0.0019  pm 

0.0021  pm 

Table  10.  Sub-range  analysis  SEPs  for  line  height  prediction  using  unknown  line 
width  of  the  0.5/1  0  W  deposition  structure  A  A.-scan  measurement  technique 
and  both  diffraction  orders  were  used  for  the  analysis. 
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4.0  OVERLAY 


Overlay  is  the  measurement  of  registration  of  two  adjacent  layers  to  each  other.  Shrinking 
design  geometries,  coupled  with  increase  in  processing  steps  required  for  the  manufacture 
of  semiconductor  devices,  are  placing  an  ever  increasing  demand  on  overlay  tolerances. 
The  overlay  budget  for  0.35  pm  structures  is  70  nm  dictating  the  precision  of  the  overlay 
measurement  instrument  to  be  7  nm  or  better.  This  requirement  decreases  to  5  nm  for  0.25 
pm  feature  sizes.  Optical  microscopy  techniques  have  been  the  dominant  measurement 
technique  to  date,  however,  it  is  not  clear  that  they  will  be  able  to  meet  the  challenge  as 
the  feature  size  approaches  0.25  pm.  The  two  major  sources  of  error  in  overlay 
measurement  are  wafer  induced  shift,  or  WIS,  due  to  overlay  mark  asymmetry,  and  tool 
induced  shift,  or  TIS,  due  to  measurement  tool  asymmetry.  It  is  clear  that  in  order  to 
achieve  overlay  measurement  precision  of  better  than  5  nm  any,  and  all,  sources  of  error 
must  be  removed. 

We  proposed  to  test  diffraction  based  techniques  for  overlay  measurement  and  present  the 
conclusion  of  our  theoretical  study.  The  structure  under  consideration  is  shown  in  Figure 
8  and  consists  of  a  developed  photoresist  grating  piaced  on  top  of  an  etched  Si  grating. 
The  region  in  the  space  of  this  grating  is  filled  with  Si02.  The  periods  of  the  two  gratings 
are  the  same  The  probe  beam  is  from  a  He-Ne  laser  operating  at  632.8  nm  and  a  2-0 
scatterometer  arrangement  is  used  for  obtaining  the  measurement.  The  diffracted  power  in 
both  the  +1  and  -1  order  beams  is  measured  as  the  angle  of  incidence  is  scanned.  The 
quantity  of  interest  here  is  K(0),  which  is  the  ratio  of  the  power  diffracted  into  the  +1 
order  at  an  angle  of  incidence  0,  to  the  power  diffracted  into  the  - 1  order  due  to  an  angle 
of  incidence  of  -0.  Because  the  period  of  the  two  gratings  is  the  same,  they  can  be 
considered  to  be  a  single  periodic  structure  and  K(0)  is  a  measure  of  symmetry  of  this 
periodic  structure 

In  the  absence  of  registration  errors,  the  photoresist  grating  will  be  patterned  exactly  on 
top  of  the  Si/Si02  grating,  and  a  symmetrical  periodic  structure  will  be  created.  This  will 
result  in  K(0)  =  1  because  the  structure  will  appear  the  same  at  an  angle  of  incidence  equal 
to  0  or  -0  Registration  errors  in  the  patterning  of  the  top  photoresist  grating  will  lead  to 
an  offset  in  the  top  grating  as  opposed  to  the  bottom  grating.  This  will  create  an 
asymmetry  in  the  periodic  structure  which  is  shown  in  Figure  8  This  will  result  in  the 
deviation  of  K(0)  from  I .  Figure  9  is  a  theoretical  plot  of  K(0)  as  a  function  of  0  for 
different  overlay  errors.  It  can  be  seen  that  the  curves  differ  significantly  even  for  a  small 
value  of  overlay  error.  Using  the  RCWT  model  we  gencrated  theoretical  plots  of  K(0)  for 
overlay  errors  ranging  from  -40  nm  to  +40  nm  with  a  step  size  of  5  nm.  These  plots  were 
input  into  a  PLS  based  calibration  model,  and  an  estimator  of  overlay  errors  was 
generated  The  cross  validated  SEP  of  calibration  was  less  than  0.5  nm,  significantly  less 
than  the  5  nm  requirement  described  previously. 
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Overlay  Photoresist 


Si/Si02 

Grating 


Figure  8:  Illustration  of  the  structure  used  to  simulate  application  of 
scatterometry  to  monitor  overlay,  consisting  of  a  Si/SiC>2  grating  with  a  grating  of 
photoresist  on  top  and  slightly  displaced  laterally. 


Figure  9:  Plots  of  K(0)  versus  angle  of  incidence  for  five  different  overlay  errors 
illustrating  the  sensitivity  of  the  scatterometry  technique  to  monitor  overlay. 
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5.0  INITIAL  INSTRUMENT  DESIGNS 


Based  on  the  results  described  above,  twc  scatterometer  configurations  appear  most 
feasible  for  applications  to  characterize  line  shapes  in  a  cluster  tool  environment  and 
overlay.  One  or  both  of  these  will  be  constructed  in  a  Phase  II  effort  as  discussed  in  the 
nexl  section. 

Figure  1 0  illustrates  an  arrangement  which  represents  a  no-moving-parts  version  of  the  \- 
scan  scatterometer  of  Figure  3.  Note  that  the  materials  involved  in  this  system  are  vacuum 
compatible,  and  this  permits  locating  the  scatterometer  in  the  load  lock  of  a  cluster  tool 
system.  A  fiber  optic  directs  white  light  to  illuminate  the  sample.  The  0-order  component 
which  is  specularly  reflected  from  the  sample  is  directed  back  through  the  same  fiber  into  a 
miniature  spectrometer  equipped  with  a  1024  element  ccd  detector  array  to  provide  rapid 
spectral  resolution  of  this  component  The  detector  arrays  shown  on  both  sides  of  the 
fiber  characterize  the  higher-order  diffraction  from  the  sample.  Different  wavelength  1  st- 
order  diffraction  components  will  illuminate  different  detectors  on  the  arrays,  thus 
providing  spectral  resolution.  This  arrangement  neglects  the  overlap  of  a  1st -order 
component  of  one  wavelength  and  a  2nd-order  component  of  half  the  same  wavelength. 
However,  we  do  not  see  this  to  be  a  problem  because  of  the  statistical  approach  we  are 
using  for  analysis  and  the  fact  that  this  "overlap"  issue  is  -epeatable.  In  addition,  the  2nd- 
order  terms  typically  are  significantly  less  ini  ;nse  than  the  1  st-order  terms  and  thus 
represent  a  minor  changes  to  the  signal  levels  in'  olved. 


To  Light  Source  and 
Miniature  Spectrometer 


Figure  10:  Arrangement  of  a  no-moving-parts  X-scan  scatterometer  which  is 
vacuum  compatible 
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The  exact  number  of  these  detectors  that  is  required  is  not  known  now,  but  we  estimate 
16  on  each  side  of  the  fiber  will  be  sufficient.  This  is  based  on  collaborative  work  with 
AT&T  in  characterizing  0.5  pm  pitch  structures.  In  this  work  we  obtained  a  complete 
data  set  of  approximately  300  points  to  determine  the  lineshape  of  the  structure  and 
validated  this  with  SEM  characterization.  We  then  performed  analysis  using  fewer  and 
fewer  data  points  until  the  predicted  lineshape  began  to  differ  from  that  predicted  with  the 
full  data  set.  We  found  that  12  data  points  provided  estimates  of  the  lineshape  that  were 
effectively  the  same  as  the  estimate  obtained  using  the  frill  data  set 

Figure  1 1  illustrates  a  scatterometer  configuration  that  is  similar  in  principle  to  that  of 
Figure  10,  but  which  is  appropriate  for  characterizing  actual  device  patterns  as  opposed  to 
line/space  test  patterns.  In  this  arrangement  the  sample  is  illuminated  with  white  light  via 
a  fiber,  and  the  specularly  reflected  0-order  component  is  analyzed  by  the  miniature 
spectrometer  as  before  The  higher-order  diffraction  from  the  device  pattern  is  directed 
through  a  lens  to  a  color  camera  or  an  imaging  spectrometer  to  provide  a  spectrally- 
resolved  measurement  of  these  components  Using  a  camera  is  simpler  and  faster,  but  it 
provides  lower  spectral  resolution.  Lower  resolution  may  not  be  an  issue  because  of  the 
statistical  analysis  used;  this  would  need  to  be  investigated.  The  lens  could  be  located 
inside  the  vacuum  load  lock  of  a  cluster  tool,  and  the  remaining  parts  of  the  system  would 
be  in  atmosphere 


Figure  1 1 :  Arrangement  of  a  no-moving-parts  X-scan  scatterometer  which  can 
characterize  device  patterns 
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Characterizing  actual  device  patterns  is  more  direct  than  characterizing  line/space  test 
patterns,  and  in  this  sense  it  is  more  appealing.  However,  it  is  not  possible  to  model  the 
two-dimensional  diffraction  from  devices,  and  this  precludes  basing  the  technique  on  a 
firm  theoretical  background.  In  particular,  training  the  statistical  analysis  technique  would 
be  based  on  an  empirical  approach  in  which  diffraction  from  a  number  of  samples  is 
characterized,  and  the  lineshape  or  overlay  of  the  samples  is  also  well  characterized.  The 
application  also  requires  the  pattern  to  have  a  high  degree  of  periodicity,  and  this  is  often 
the  case  in  advanced  devices  such  as  microprocessors  and  memory  arrays.  This  approach 
has  been  used  at  the  University  of  New  Mexico  to  accurately  characterize  trench  depths  of 
16  MB  DRAM  structures  and  planarization  of  SRAM  device  patterns. 


I 

I 
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6.0  CONCLUSIONS 


We  have  shown  that  two  scatterometry  techniques  can  be  applied  to  measure  the  line 
shapes  of  etched  poly-Si  and  deposited  W  structures  having  0.20  pm,  0.35  pm  and  0.50 
pm  CDs  with  errors  that  are  well  within  acceptable  tolerance  levels.  These  results  can  be 
improved  using  the  sub-range  technique  described  In  addition,  overlay  can  be  measured 
with  less  than  0.5  nm  error. 

Simulated  scatterometer  data  was  calculated  using  rigorous  coupled  wave  diffraction 
theory.  Analysis  was  performed  using  statistical  techniques  (PLS)  to  provide  a  simple 
method  to  determine  the  performance  that  can  be  expected  from  the  scatterometer 
techniques.  In  addition,  these  simulation  techniques  represent  a  significant  simplification 
compared  to  generating  samples  and  characterizing  them;  this  will  be  performed  in  a  Phase 
II  effort.  We  have  used  these  same  simulation  techniques  in  similar  applications  in  which 
experimental  results  were  in  close  agreement  with  simulated  results.  These  other  results 
involved  CDs  between  0. 12  pm  to  2.0  pm. 

This  work  has  laid  the  foundation  for  developing  more  sophisticated  scatterometer 
arrangements  by  showing  the  feasibility  of  the  two  basic  techniques,  2-0  and  L-scan,  to 
measure  lineshape  and  overlay  Using  thi'  background,  more  advanced  systems  can  be 
designed  and  constructed  to  address  specific  issues,  such  as  cluster  tool  applications 


7.0  PHASE  II  PLAN 


We  plan  to  propose  a  Phase  II  effort  which  will  involve  constructing  scatterometers  and 
implementing  them  on  an  etch  or  deposition  reactor  to  pursue  in-sitv  measurements. 
Presently  we  are  having  discussions  with  Sandia  National  Laboratories  (SNL)  and  Lam 
Research  in  this  regard  SNL  is  presently  establishing  the  capability  to  process  0.5  pm 
geometries  However,  smaller  geometries  would  require  ex-situ  measurements  using 
samples  which  would  be  obtained  elsewhere.  One  possibility  we  are  exploring  is  using  e- 
beam  written  patterns  in  photoresist  to  fabricate  smaller  geometries  with  ion  beam  milling. 

We  will  also  involve  the  University  of  New  Mexico  as  part  of  our  Phase  II  effort  because 
of  their  background  in  demonstrating  scatterometer  applications  in  microelectronics 
manufacturing.  They  will  provide  the  modeling  required  for  the  fundamental 
understanding  of  processes  as  well  as  the  optical  expertise  required. 
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Scatterometry  Applied  to 
Microelectronics  Processing 

J.R.  McNeil,  S.S.H.  Naq  vi,  SAL  Caspar,  K.C.  Hickman,  K-P.  Bishop,  LJVL  Milner,  RJL  Knikar,  G .A  Petersen 

Center  for  High  Technology  Materials,  University  of  New  Mexico,  Albuquerque,  New  Mexico 


Scatterometry  is  a  promising  new  optical  metrology 
method  for  nondestructive  rapid  evaluation  of  many 
physical  quantities  of  critical  interest  in  microelectronics 
fabrication.  Among  these  are  latent  image  CDs  in 
undeveloped  resist,  lithography  tool  focus  and  exposure 
parameters,  etched  structure  CDs,  photomask  and  phase 
shift  mask  characteristics,  and  metal  grain  sizes.  This 
paper  provides  a  general  review  of  scatterometry,  describes 
both  Instrumentation  and  methods,  and  introduces 
subjects  to  be  discussed  later  In  this  series. 


c 

tePritical  dimensions  (CDs)  of  0.8  pm  are  common  for  high 
performance  microelectronic  devices.  Advanced  fabrication 
processes  already  involve  CDs  of  0.5  pm  while  technology  for 
CDs  of  0.35  pm  and  smaller  is  now  being  developed  and  the 
0.1-  to  0.2-pm  range  will  be  broached  by  early  in  the  next 
decade.  This  trend,  coupled  with  increasing  automation  and 
larger  die  and  wafer  sizes,  only  hints  at  the  unprecedented  de¬ 
mands  now  being  made  in  the  area  of  fabrication  process  con¬ 
trol.  With  current  CD  error  budgets  typically  10%  to  20%  of 
the  total  CD,  the  absolute  dimensional  control  required  in  fab¬ 
ricating  advanced  devices  is  0.05  pm  or  smaller. 

Unprecedented  Process  Control  Demands  an 
Unprecedented  Metrology 

Ideally,  a  metrology  technique  should  be  noncontact,  nonde¬ 
structive.  rapid,  quantitative,  simple,  and  applicable  in  situ.  An 
attractive  alternative  metrology  technique  that  fits  we  '  with 
these  objectives  has  appeared  relatively  recently.  Called  scat¬ 
terometry,  it  is  the  angle-resolved  characterization  of  light  scat¬ 


tered  from  a  surface.  The  method  may  be  used  for  characteriz¬ 
ing  structure  and  topology  as  well  as  for  controlling  fabrication 
processes.  It  is  noncontact,  nondestructive,  rapid,  and  quantita¬ 
tive.  and  in  many  cases  can  be  used  in  situ.  Scatterometry  is 
useful  in  many  areas,  including  microelectronics  and  the  pro¬ 
cessing  of  optical  elements  and  storage  media.  Its  use  can  elimi¬ 
nate  the  need  for  optical  microscope  examination  in  many  cases 
and,  unlike  microscopy,  it  is  readily  amenable  to  automated 
processing.  In  this  article,  we  will  review  scatterometry  applica¬ 
tions  that  include  characterization  of  photoresist  latent  images 
for  exposure  and  focus  control,  etched  structures,  photomasks 
and  phase  shift  masks,  and  Al-Si  grain  size. 

Basic  Principles 

Scatterometry  is  a  light  scattering  technique.  Light  scattering 
measurements  have  been  used  since  the  mid-1970s  to  charac¬ 
terize  optical  elements  [1-4]  in  terms  of  the  root  mean  square 
(rms)  roughness.  More  recently,  scatterometry  has  been  ap¬ 
plied  to  materials  and  processes  used  in  microelectronics  fabri¬ 
cation  [6-9].  The  method  can  characterize  structures  with 
dimensions  much  smaller  than  the  wavelength  X  of  the  light 
used,  in  contrast,  optical  microscopy  has  a  spatial  resolution 
comparable  to  X. 

Samples  characterized  by  scatterometry  have  a  periodic 
structure.  This  results  from  intentionally  patterning  the  sample, 
or,  alternatively,  by  considering  a  nominally  smooth  sample  to 
be  a  superposition  of  many  periodic  structures  as  described 
below.  Light  which  illuminates  such  a  sample  is  diffracted  at 
angles  prescribed  by  the  glaring  equation 

sin  ©j  +  sin  ©f  =  n XI d  ( 1 ) 

where  ©j  and  @t  are  the  incident  and  scattered  angles, 
respectively,  n  is  an  integer,  and  d  is  the  period  or  pitch  of  the 
structure.  With  scatterometry,  unlike  imaging,  the  diffracted 
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txppr  ital  arrangement  of  a  scatterometer 


Srotterometer  arrangement  for  characterizing  two-dimensional 
srnltni  patterns 


lip. lit  is  nut  combined  or  imaged  prim  m  analysis  Instead, 
inlomialimi  about  the  lineshape  ol  a  sample  is  derived  from  the 
relatiu'  intensities  nl  the  dilferent  diffraction  orders  from 
la |  1 1 ).  zero-order  and  Ist-order  diffraction  can  occur  tor  sam¬ 
ple'-  hiving  a  pitch  slightly  larger  than  A/2.  Thtis.  while  the 
iiinnhei  ol  orders  mas  be  Imnled  In  A,  the  lineshape  mlorma 
lion  is  limited  mils  In  the  sensilivits  with  which  small  sarin 
lions  in  the  intensities  ol  these  orders  can  be  detected  The 
gre.itei  the  number  ol  orders,  die  greater  the  information  avail 
able  hut.  as  sse  shall  see.  much  intorinaiimi  can  he  gleaned 
trom  the  zero  order  alone  Cnmpuler  simulations  suggest  that 
seal! ermnet i  s  e.m  be  used  in  pi ae lice  lo  characlenze  (  I  )s  as 
small  as  Z./20  Periodic  test  patterns  are  not  necessarily  le 
tpnred  ns  mnns  microelectronic  desiees  bust  pciiodic  sirue 
lines  ss Inch  dillinel  light 

Experimental  Arrangement 

I  mine  I  shosssnbn-.n  scnllci omelei  urrnnremeni  flic  luthl 


A  "2(-l  ScaMeromeler"  arrangement  especially  suited  lor 
characterizing  shoo  pilch,  small  CD  structures 


source  is  usually  a  laser  to  simplify  alignment  and  provide  ample 
signal  Typically,  the  plane  of  polarization  is  perpendicular  to  the 
plane  of  incidence  is- polarized  light )  as  this  avoids  sin taee  plas¬ 
ma  wave  coupling  m  examining  conducting  samples.  Often  the 
laser  output  is  spatially  tillered  to  provide  a  well-defined  spo!  at 
the  sample.  This  is  critical  for  allowing  measurements  close  to 
the  specularly  reflected  or  directly  transmitted  beams.  In  addi¬ 
tion.  the  polarization  of  the  detected  lighl  can  be  analyzed  using 
standard  inothinls  described  in  the  literature  [3-bj. 

Two  dimensional  scatter  patterns  are  easils  viewed  using  the 
arrangement  oflig.  2.  Here  the  sample  is  placed  al  the  center 
ot  curvature  ol  a  frosted  dome  and  is  illuimnaled  so  that  the 
scatter  pattern  is  displayed  on  the  inside  surface  ol  the  dome.  A 
camera  on  the  opposite  side  ot  the  dome  eaptuies  the  image,  a 
small  computer  is  used  tor  image  processing  and  analssis. 
there  are  no  moving  parts,  and  measurements  ate  iapnl.  l-.adi 
spot  ot  the  pattern  represent-,  a  diffraction  order  and  contains 
information  about  the  diffracting  stnictute  figure  t  illustrates 
another  arrangement,  which  is  well  suited  tor  characterizing 
short  pitch,  small  Cl)  structures  Ibis  20  Nc.itteiometei "  luo 
previously  been  employed  hv  /.nidi,  el  al.  |  101  tor  grating  chat 
actenzation.  lis  use  involves  rotating  the  sample  to  continuous 
ly  vary  the  angle  ot  incidence.  The  detector  rotates  .it  a  velocity 
such  that  the  intensity  ot  a  specific  diltrnction  oidei  is  column 
ously  monitored,  thus  providing  a  huge  daia  set  lor  that  ordei 

A  Theoretical  Overview  of  Scatterometry 

Diffraction  From  Periodically  Patterned  Surfaces 

When  .in  electromagnetic  wave  is  hr  idem  on  a  surface,  the 
light  which  i-  sc.Ulcicd  is  delcmnned  bv  the  i  ompositioii.il 
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Diffraction  characteristics  for  two  periodic  structures  having  the 
same  pitch  but  different  lineshape 


properties  <e  g.,  refractive  t rules  of  the  different  regions)  as 
well  as  the  topographical  characteristics  (e.g.,  pilch,  depth,  and 
sidewall)  of  the  surface  Properties  of  the  incident  beam  (in¬ 
cluding  wavelength,  angle  of  incidence,  and  polarization)  also 
determine  the  diffraction  characteristics  of  die  sample  Modi 
the  reflected  and  transmitted  electromagnetic  fields  carry  infor 
malum  about  the  scattering  surfau.  Numerous  metrology  tech 
tuques  te  g,,  ellipsometry.  optical  microscopy,  etc.)  exploit  this 
property  ol  scattered  fields  to  obtain  compositional  and/or 
topographical  mlnrmalion  about  the  scattering  surface.  When  a 
periodic  surface  is  illuminated,  die  scattered  fields  consist  of 
distinct  diffraction  orders  propagating  in  directions  specified 
fry  the  grating  equation  This  is  shown  schematically  m  f  ig  4, 
which  illustrates  the  diffraction  order  intensities  of  two  .true 
lures  tli.it  have  the  same  pitch  hut  different  linew  tilths 

1  lie  topographical  dependence  ol  diffraction  is  further  illus 
(rated  m  lig.  5,  which  shows  the  calculated  I st-order  dittrae 
lion  intensity  versus  linewulth  lor  a  1.5  (mi  pitch  periodic 
photoresist  structure.  The  sample  was  a  1.0  gun  photoresist 
layer  on  150  nm  ol  Sit ),  on  a  Si  substrate  The  OVI  um  beam 
was  incident  at  40  and  was  Tli  polarized.  Diffraction  charac 
tensin',  plotted  in  the  ligiue  were  calculated  by  direct  solution 
ol  Maxwell's  equations  I  I  I  I  Details  will  fie  given  in  the  last 
article  ol  itus  scaltcroinctrv  'cues 

Topology  of  "Smooth"  Surfaces 

The  relation  between  electromagnetic  scattering  anil 
smooth  surface  topology  has  been  studied  for  manv  years, 
oiigmallv  in  connection  with  radar.  In  the  case  of  a  clean, 
perfectly  reflecting  surface  where  height  variations  are  much 
smallei  than  the  wavelength,  the  situation  is  simplified  I -'of 
lowing'  ( 'lunch's  treatment  |  1 ,2  |  the  sample  is  modeled  as  an 
mtiinic  niimbei  of  surfaces,  each  having;  structure  will)  one 
spatial  frequency  and  dittrae  ling  light  at  one  angle  ui  cording 
to  l\q  (  li  The  sm face  is  assumed  to  be  sufficiently  smooth 
such  that  onls  1st  order  diffraction  is  significant 


Model  dotn  dlustrotes  the  dependence  of  1  st-order  diffracted  power 
on  the  imewidlh  ol  periodic  structure  having  a  I  5-pm  pitch 
Conditions  simulated  were  Si  substrate  with  350  nm  of  SiO?,  1  pm 
photoresist,  ?.  ■  633  nm,  beam  TE  polarized  and  incident  at  30 


Vector  scattering  theories  desc  ribe  the  differential  light  seat 
tcrdl  as 

s 

Idl  C 

— - -  -  Q((-> ,<P|, I’< p.c| >  ( 2 ) 

I,  citi),  X* 

where  C  is  a  constant.  is  the  intensity  of  the  incident  light, 
and  dm,  is  the  solid  angle  of  the  detection  system.  The  quantity 
U  m  Hq.  (2),  called  the  “optical  factor",  is  independent  ol  the 
surface  structure  and  is  a  function  of  the  angles  of  incidence 
(W,  0,).  si  alter  angles  ((-)  ,<J)J,  complex  index  of  retraction  N  ol 
the  surface,  and  polarization  states  of  the  incident  and  scattered 
light,  x,  und  X,.  respectively.  T  he  “surface  factor"  I’tp.qi  is  the 
power  spectral  density  <PSP|  of  the  surface  roughness  I  bis  is 
the  output  of  the  scatterometer  measurement  and  is  the  tunc 
turn  that  describes  the  surface  strueture 

If  the  surface  is  m  the  ,V  T  plane,  and  z't.Y,  )T  is  the  surface 
height  variation  relative  to  that  plane,  the  I’SD  is  given  by 


I’(  Pdf) 


\\dxdy  ,.*<!«  "f»>Z(Xly) 


(T) 


where  A  is  the  area  illuminated,  and  p  and  q  are  [he  surface 
spatial  frequencies  ui  the  ,V  and  Y directions,  respectively  In 
other  words,  the  PSI)  is  the  average  squared  magnitude  <d  the 
two -dimensional  I-ouricr  transform  of  the  surface  roughness 
It  is  convenient  todesciibc  the  surface  roughness  in  terms  of 
the  rzm  roughness  a,  which  is  given  in  terms  of  the  instrument 
bandwidth  and  modulation  transfer  function  (M  ITT.  Mlp.ql  as 
Different  values  of  o  will  result  if  the  integral  limits  ti.e  , 
p.~.  >1*.. 

<T‘  -  Jd[)  Jdq  M(p.q)Plp.q)  (-li 

l>™„  >i™„ 

bandwidth  I  or  M  IT  in  the  integral  change  II  rm s  roughness 
values  from  examining  die  same  surface  using  diffeient  msini 
incuts  .lie  to  Ire  compared,  tfie  haudw tilths  and  M  IT's  should  he 
file  same  |4|  In  the  ease  ol  an  isotropic  surface  sp  tm,  p  utul 
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Figure  6  _ 


Figure  7 
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Changes  m  developed  photoresist  CD  resulting  from  variations  in 
underlying  poly  Si  film  thickness  The  scatteiometercontrolled 
I'xposme  dose  provided  nearly  a  factor  of  four  improvement  in  CD 
control  compared  to  using  fixed  exposure  time  The  linewidth  and 
pound  were  nominally  I  pm  and  ?  pm.  respectively 


q  ;in*  iduntKal.  .uul  otil\  one  of  the  spatial  wau-k-nptfis  is  need 
oil  to  describe  the  lateral  dimension  of  sample  slim  tun*. 

Applications  Overview 

A  Il'Ml'W  (it  sC.ltteiOIlk’IIV  .1|’|>  I U atlolls  I1*  pil’SCIdeil  lietc 
Tuttlier  details  will  appear  in  subsequent  articles  ol  this  senes 
lot  applications  m  micioelectrnmcs.  sample  surfaces  ■  an  he  a 
lest  pallet  n  m  pan  ul  l  lie  device  hemp  lain  lea  led  I  lie  i  llu  mi 
ikk i nr  spot  ean  he  as  small  as  til  pm  i  l/e'  diamelei  i.  so  llial  a 
minimum  ul  ten  m  so  lines  ul  the  peinulie  sluieuite  aie  illiimi 
iiali’il  | (> |  Snell  a  small  lest  palletn  ean  he  loealed  in  llu*  scribe 
lines  o|  tin-  walci 

Photoresist  Latent  Image  Characterization 

I  n  1 1 1  In  >l:  i  apliv .  e  epos  hi  e  and  I  oeus  are  hiu  moot  the  van 
allies  llial  control  the  shape  ot  the  developed  pholoicsisi  pal 
let  it.  flak  ini?  subsequent  to  exposure  also  mlltienees  the 
i esu 1 1 in o  shape  ol  the  deudoped  resist  and  thus  represents  an 
in lil i ( : o n .il  vat  ilile  hifli  must  lie  accounted  lot  m  process 
evaluation  It  .  prelerahle  to  decouple  these  elleets  and  Hide 
pcndcnllv  de  t mute  the  appiopnatr  process  w  itulou  lor  each 
ol  the  evpoMiie.  hake,  and  develop  steps.  Scallcrniueti  v  ean  he 
used  elleettveh  in  this  rcgaiil  to  momloi  the  latent  mi. me  pitot 
lo  hake  and  develop 

In  photoresist,  a  latent  linage  is  loi tiled  as  the  photo. telive 
compound  tl’Al'i  concentration  changes  upon  exposme.  A 
'.’lalum  mi  the  photomask  can  gem  t.ile  a  resist  pattern  with  a 
pet  lodii  I’  At  A  mu  cut  i  at  ion  pi  ol  tie  that  ill  1 1  r.u  Is  he  lit  vv  hen  d 
liimmaled  |l.'.l  *|  llu-  dilhavtion  chat  aetet  is|  u  s  o|  latent 
lilt. tee  I'l.itine  .  ,  an  he  iipnmuslv  modeled  I  list.  Ittlioei.iphv 
simitlatimi  tools  such  as  SAMI’I  I  |  I  1 1  and  l’l«  )l  I  I  1 1  [Is!  atc 


Variation  of  1  sto'der  diffraction  intensity  from  o  photoresist  Intent 
imoqe  for  different  stepper  locus  conditions 


used  to  predict  the  PAC  eoneettiration  inside  (he  photoresist  as 
a  Inaction  ol  exposure  conditions  ulose,  loeusi  and  substrate 
and  resist  parameters  (thickness,  composition).  Then.  rigorous 
coupled  wave  theory  |  l(i|  is  employed  to  predict  the  amount  ol 
incident  probe  beam  power  diffracted  into  the  various  orders. 

A  significant  problem  occurs  in  optical  lithographs  because 
the  exposure  dose  is  a  function  ot  the  underlying  lilm  structure. 
Since  the  underlying  films  can  vary  Irom  lot  lo  lot.  water  to 
waler.  ot  even  die  to  die  on  a  wafer,  the  use  ot  fixed  exposure 
time  tdoset  the  piesent  techmi|iie  used  lor  CD  control 
ean  adversely  alfeet  critical  dimension  control.  Hv  using  seal 
(eroinetry  to  momloi  the*  latcni  image.  CD  sanations  can  he  te 
duccil  | h.d. I  2 1  by  nearly  lour  tunes  d  ig  fi).  These  results  arc 
applicable  to  any  transparent  lilm  material;  poly  St  was  used  in 
llus  work  lo  provide  the  electrical  measurements  necessary  lot 
CD  verification. 

Diffraction  Irom  latent  image  gratings  can  also  he  used  to 
determine  optimal  focus  conditions.  Theoretical  simulations 
have  shown  | ]  that,  lot  alls  exposure  dose,  the  power  dll 
track’ll  mlo  the  1st  order  is  at  a  peak  for  optimal  locus  cotkli 
nous,  regardless  ot  the  composition  ol  layers  underneath  the 
photoresist.  Tins  was  lirst  demonstrated  experimentally  hy 
Adams  [  I  (|  Milner  et  id.  |‘l|  established  that,  tortile  case  ol 
photoresist  on  crystalline  Si.  the  1st  and  2nd  urdei  ditlruc 
Mon  intensities  reach  maxima  at  optimal  locus  However,  te 
cent  experimental  results  Tor  photoresist  on  u  Si  ami  polv  Si 
show  that,  while  the  1st  older  intensity  always  peaks  at  upn 
mal  locus,  i he  2nd  order  intensity  is  dependent  upon  the  im 
deriving  layers  and  mav  or  may  not  peak  at  optimal  locus 
Thus,  it  tv  important  to  distinguish  between  the  various  dil 
traded  orders  and  detect  each  ot  them  separately 

liguie  7  illustrates  experimental  and  theoretical  1st  ordet 
intensity  variation  versus  locus  Im  photorcMsi  on  </  Si  The 
zero  delocus  position  is  at  the  top  ol  the  photoresist  The  dil 
terenee  between  the  peak  location  of  the  two  curves  is  believed 
lo  he  due  to  an  ertot  m  the  ability  ol  the  stepper  to  locate  the 
top  id  the  pholoicsisi  Niihseipi-  ut  STM  measnieitienis  id  the 
developed  photoresist  matings  unhealed  that  the  opum.  1  Im  us 
location  was  at  II.  I  |iiu.  the  same  as  the  peak  ul  the  espei  mien 
tal  v.  lit  v i’ m  the  ligute 
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Anolysis  strategy  used  to  predict  sample  etch  depth,  theoretical 
scatterameter  data  were  used  to  tram  statistical  techniques  and 
experimental  scatteromcter  data  were  subsequently  input  to  predict 
the  depth  of  a  sample 


Scatterometry  for  Etched  Structure  Characterization 

Scatterometry  cart  he  used  to  analyze  etched  structures 
rapidly  and  nondestructively  ( 1 7  IX|  One  example  of  scat- 
tcmineter  characteri/.atmn  ol  the  depth  id  an  etched  structure 
involved  patterns  that  had  12-pm  pitch  and  I  l-pm  wide  lines 
etched  in  crystalline  Si.  Due  to  the  large  pitch,  the  diffraction 
patterns  from  these  structures  contained  many  orders,  and  the 
intensity  of  each  order  was  sensitive  to  the  depth  of  the  lines. 
Statistical  techniques  were  used  to  analyze  the  patterns,  as  il¬ 
lustrated  schematically  in  big  H.  The  goal  was  to  predict  the 
etch  depth  from  the  scatter  characteristics.  I’rocedtirally,  this 
involved  the  use  ol  a  neural  network  and  a  partial  least  squares 
method  1 11)  2  1 1  hirst.  the  coupled  wave  model  generated  a  set 
ol  theoretical  data  that  was  used  tor  calibration  and  for  "tram 
mg'"  the  statistical  techniques  Then,  experimental  data  from 
the  scatter mietcr  measurements  was  input  to  the  statistical  op 
(•rations  to  predict  the  sample  depth  Results  were  very  good, 
l  ot  lourteen  samples,  each  nominally  '>00  run  in  depth,  the 
neural  net  prediction  had  an  average  error  ol  XX7  *  4  <Ki7 .  or 
XI)  nm  *  17  mn  The  partial  least  squares  prediction  tightened 
this  to  an  average  error  ol  O.X7r  *  0.17,  or  7  run  i  1  nm  Sam 
pie  depths  were  verified  using  a  scanning  force  microscope 
i  STM  I.  which  had  an  error  ol  approximately  17  It  was  not 
necessary  to  fabricate  and  characterize  calibration  samples 
The  calibration  data  that  were  generated  from  the  model  arc 
based  on  Maxwell’s  equations  and  rigorous  dillraction  theory 
and  thus  can  lie  considered  absolutely  accurate.  Hus  same  ap 
proach  ctm  be  extended  to  predict  sidewall  characteristics 

Chrome-on-Glass  Photomask  Grating  Metrology 

Act  urate  luiewidlh  metrology  ol  photomasks  is  essential  for 
production  o|  devices  with  I-  it urc  sizes  less  than  0  s  pm  A 
NISI  <  I !  S  National  Institute  ol  Standards  and  Technology! 
Inn  width  standard  now  exists  lor  O')  pm  lines  ol  t  hromc  on 
glass  However,  linewidth  standards  are  nonexistent  lor  smaller 
f  "IK  ol  i  hr oine  on  glass  or  am  other  material,  lor  that  matter 


Model  and  experimental  dola  loi  the  zero-order  diffraction  intensity 
from  a  line  ctructuie  with  a  ICT|irr  pitch,  5)im  width,  and  0  371  | i/n 
depth  measured  using  the  "2H  hcotteromelor  "  arrangement 


Scatterometrv  provides  a  simple  and  accurate  method  for 
determining  line-widths  of  photomask  gratings  In  this  ap 
plication,  an  incident  laser  beam  produces  a  scattered  field 
consisting  of  multiple  reflected  and  transmitted  orders.  A 
simplified  scalar  model  of  diffraction  indicates  that  the 
amount  ol  power  diffracted  into  the  ttansmittcd  zero  order 
should  he  inversely  proportional  to  the  Imewtdth  The 
linewidth  ol  the  photomask  grating  may  he  predicted, 
therefore,  from  a  measurement  of  intensity  ol  the  transmit 
led  zero  order  |K|.  To  rigorously  establish  this  inverse  rela 
tionship,  the  chrome  ori  gins-  grating  may  be  modeled  as 
an  XI)  nm  thick  chrome  la'  et  deposited  on  a  glass  -ubsttale 
with  an  I  X  tun  thick  antircflectmg  layer  ol  chrome  oxide 
Application  ol  rigorous  coupled  wave  theory  |lh|  con 
firmed  a  monotonic  relationship  between  the  linewidth  ol 
the  grating  and  the  transmitted  zero  order  pawn.  A 
"  lookup  table "  can  be  used  to  give  the  linew  idth  lor  a  mea 
sured  zero-order  intensity, 

Hxperiments  verified  that  the  linewidth  obtained  with  scat 
terometrv  was  cotnparahle  to  that  obtained  using  other  opt  ten! 
metrology  instruments  One  mask  having  six  dilfereni  siibmi 
cron  linewidth  gratings  was  sent  to  several  companies  to  delet 
mine  instrumciu-to-uistruinent  variation  in  CD  measurements 
No  consensus  was  obtained  lor  the  linewidth  values  using  the 
different  instruments.  Scatterometry  piovided  Cl)  measure 
incuts  that  were  comparable  to  those  obtained  using  the  van 
nos  metrology  instruments  In  some  cases,  the  results  differed 
hv  less  than  17  (e.g.,  results  Iron)  a  Namibia'  1 1  111 )  The  long 
term  t  I  X  mus  l  repeatability  ol  the  measurements  was  lottml  to 
lx-  approxunatelv  'M>7 

Line  Edge  Roughness 

Mask  patterns  can  have  both  random  stnu  lure  and  periodi; 
cnoi  on  line  edge  .  The  lormer  pioduces  an  isotropic .  dillu  -- 
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Power  spectral  density  (PSD)  characteristics  of  Al-Si  material  deposited 
on  Si  wafers  at  diffemt  c  eposition  temperahtres. 


scatter  und  the  later  produces  diffraction,  both  of  which  are 
superimposed  on  the  diffracted  pattern  of  the  desired  line 
structure.  These  two  undesirable  components  of  scattered 
light  can  be  treated  separately,  but  preliminary  work  has 
focused  on  a  differential  measurement  technique  that  in¬ 
cludes  both. 

The  dome  technique  of  Fig.  2  was  used  to  measure  the  two- 
dimensional  scatter  from  mask  patterns  that  nominally  have  a 
one-dimensional  periodic  structure.  For  example,  perfect  struc¬ 
ture  on  the  mask  would  diffract  through  a  diameter  of  the  dome 
at  several  points  determined  by  the  number  of  diffraction  or¬ 
ders  that  occur.  For  each  measurement,  the  data  are  transferred 
to  a  computer  and  image  processing  techniques  subtract  the 
ideal  diffraction  pattern  from  the  experimentally  obtained  pat¬ 
tern.  The  contribution  of  the  residual  scatter  is  integrated  and 
normalized  to  the  integral  of  the  total  scatter  in  the  experimen¬ 
tal  pattern.  This  defines  a  figure  of  merit  (FOM)  as 


FOM  = 


f  Scattered  light  exclusive  of  ideal  diffraction 
J  All  scattered  light 


Three  photomasks  were  examined:  one  which  had  good 
quality  lines,  one  that  had  a  sinusoidal  edge  structure  intention¬ 
ally  superimposed  on  the  line  edges,  and  one  that  was  inten¬ 
tionally  etched  to  produce  significant  line  edge  roughness.  The 
resulting  FOM  values  for  the  masks  were  3.57,  3.03,  and  2.56, 
respectively,  suggesting  that  scatterometry  is  a  promising 
quantitative  technique  for  mask  characterization. 

Phasa  Shift  Masks 

Phase  shift  photomasks  have  recently  attracted  a  great  deal  of 
attention.  Selective  changes  in  the  phase  of  the  light  as  it  tra¬ 
verses  the  phase  mask  improves  the  resolution  of  the  imaged 
lines  which  are  printed.  Accurate  measurement  of  the  phase 
shifts  resulting  from  thickness  variations  is  presently  a  prob¬ 
lem.  In  line  with  this,  mask  structures  were  examined  using 
the  “20  Scatterometer."  The  results  are  illur,,uted  in  Fig.  9 


Experimental  and  model  results  are  plotted  for  the  zero-order 
intensity  diffracted  from  a  10-prn  pitch.  5-pm  linewidth 
structure,  which  was  ion  beam  milled  to  a  depth  of  approxi¬ 
mately  0.37  pm.  The  angle  of  incidence  was  varied  over  ± 
50°.  The  dip  in  the  solid  curve  at  0°  results  from  the  detector 
blocking  the  incident  beam.  The  difference  between  the  ex¬ 
perimental  and  modeled  curves  is  2.3%.  Modeled  data  for 
other  depths  indicate  that  the  variation  of  the  zero-order  light 
intensity  is  quite  sensitive  to  the  etched  depth. 

Nominally  "Smooth"$urfocai 

Figure  10  illustrates  the  power  spectral  density  characteris¬ 
tics  for  five  samples  (courtesy  of  Intel)  of  Al-Si  alloy  materi¬ 
al,  which  were  deposited  at  different  substrate  temperatures 
between  ambient  and  approximately  350°C.  The  area  under 
the  PSD  curves  increases  and  the  secondary  peaks  shift  to 
lower  spatial  frequencies  (higher  spatial  wavelengths)  as  tem¬ 
perature  increases.  The  area  increase  is  indicative  of  increas¬ 
ing  sample  roughness,  as  expected,  and  the  peak  shift  is 
indicative  of  increased  grain  size  of  the  material.  Using  SEM 
techniques  to  determine  grain  size,  consistent  results  were 
obtained  when  compared  to  rms  roughness  (98%  correlation), 
although  the  SEM  characterization  required  the  samples  to  be 
broken  and  etched.  Similar  results  have  been  observed  for  Al- 
Cu  and  Al-Si  (from  other  manufacturers)  deposited  at  differ¬ 
ent  substrate  temperatures,  levels  of  substrate  bias,  and  with 
antireflection  (AR)  layers.  Samples  of  Al  have  also  been 
characterized  using  gTating  coupling  to  surface  plasma  waves 
to  determine  the  electrical  properties  of  the  films  [22]  This 
application  of  scatterometry  is  compatible  with  the  industry 
trend  to  build  characterization  directly  into  sputter  systems  or 
cluster  tool  arrangements. 

In  a  similar  manner,  we  have  characterized  other  materials 
including  CVD  W  and  silicides.  For  example,  the  effect  of 
oxide  growth  on  the  surface  properties  of  WSi3  is  easily  detect¬ 
ed  using  scatterometry.  The  as-deposited  material  is  smooth, 
and  the  PSD  peaks  at  small  spatial  wavelengths,  outside  the 
measurement  bandwidth  of  the  scatterometer.  Subsequent  oxi¬ 
dation  of  the  surface  at  elevated  temperatures  promoted  grain 
growth  in  the  material.  The  rms  roughness  of  the  sample  in¬ 
creased  by  a  factor  of  four,  and  the  sample  PSD  plot  peaked  at 
large  spatial  wavelengths  (5  pm)  This  is  consistent  with  results 
from  microscopic  examination  of  the  samples. 

Conclusions  and  Forecast  of  Future 
Applications 

The  scatterometry  applications  reviewed  above  illustrate 
the  technique’s  wide  range  of  applicability  in  providing  a 
nondestructive,  rapid,  quantitative  process  monitor.  Although 
the  applications  discussed  range  from  lithography  control  to 
etched  structure  monitoring  to  metal  grain  size  determination, 
these  only  begin  to  utilize  the  technique's  capability.  Mea¬ 
surements  can  also  be  performed  using  multiple  wavelengths 
and  polarizations  to  yield  even  greater  information.  As  a  di¬ 
mensional  metrology  tool,  scatterometry  will  meet  industry 
requirements  well  into  the  future.  Visible  light  enables  scat¬ 
terometry  characterization  of  structures  that  have  0.15-pm 
CDs  and  smaller;  use  of  ultraviolet  sources  can  extend  these 
limits  even  further. 

Scatteromcters  ol  the  form  shown  in  Fig.  1  are  commercially 
available  from  TMA  Technologies,  Bozeman.  MT.  In  addit.on. 
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